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Abstract 

We consider the non-relativistic four-boson system with short-range forces and large scattering 
length in an effective quantum mechanics approach. We construct the effective interaction potential 
at leading order in the large scattering length and compute the four-body binding energies using the 
Yakubovsky equations. Cutoff independence of the four-body binding energies does not require 
the introduction of a four-body force. This suggests that two- and three-body interactions are 
sufficient to renormalize the four-body system. We apply the equations to 4 He atoms and calculate 
the binding energy of the 4 He tetramer. We observe a correlation between the trimer and tetramer 
binding energies similar to the Tjon line in nuclear physics. Over the range of binding energies 
relevant to He atoms, the correlation is approximately linear. 
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I. INTRODUCTION 



Effective theories are ideally suited to describe the low-energy properties of physical sys- 
tems in a model-independent way. They can be applied to any system that has a separation 
of scales, which can be a fundamental property of the underlying theory or simply a kine- 
matical suppression. The long-distance degrees of freedom must be included dynamically 
in the effective theory, while short-distance physics enters only through the values of a few 
coupling constants, often called low-energy constants. Effective theories are widely used in 
many areas of physics. Recently, a considerable effort was devoted to applying effective field 
theories in nuclear and atomic physics. For overviews of these programs, see e.g. Refs. 0, 0| 
and 0, 0| . If there is no exchange of massless particles, any interaction will appear short 
ranged at sufficiently low energy. One can then use a very general effective theory with 
short-range interactions only to describe the universal low-energy properties of the system. 
Such a theory can be applied to a wide range of systems from nuclear and particle physics 
to atomic and molecular physics. 

Most previous work in this area was done using a non-relativistic effective field theory 
(EFT) with contact interactions. Particularly interesting are few-boson systems with large 
scattering length. They are characterized by an unnaturally large two-body scattering length 
a which is much larger than the typical low-energy length scale I given by the range of the 
interaction. Such systems display a number of interesting effects and universal properties 
that are independent of the details of the interaction at short distances of order I. If a > 0, 
e.g., there is a shallow two-body bound state with binding energy B 2 = h 2 /(Ma 2 ) + 0(l/a), 
where M is the mass of the particles. Low-energy observables can generally be described 
in a controlled expansion in l/\a\. In the two-boson system, the effective theory reproduces 
the effective range expansion (cf. Refs. 0, 0) but the structure of the three-boson system 
with large scattering length is richer. In Refs. 0, H| , it was found that both two- and three- 
body contact interactions are required at leading order for the consistent renormalization 
of the three-body system. Interestingly, the renormalization group behavior of the three- 
body interaction is governed by an ultraviolet limit cycle. This implies that at leading 
order in l/\a\, the properties of the three-boson system with large scattering length are not 
determined by two-body data alone and one piece of three-body information (such as a three- 
body binding energy) is required as well. In the EFT, this information can conveniently be 
parameterized by the three-body parameter A* introduced in Refs. 0, Hf. These general 
findings confirm and extend previous work by Efimov who derived many general features of 
the three-body problem with large scattering length [il. Eof. 

While the EFT formulation has been very successful, it is not the only possible for- 
mulation of an effective theory for this problem. Lepage has advocated the framework of 



non-relativistic quantum mechanics with an effective interaction potential The contact 
operators in the field theory are the replaced by an "effective potential" built from smeared 
out 5-function potentials and derivatives thereof. In the case a > 0, this approach has been 
applied to the three-boson system by Wilson l3| and Mohr 0] . They confirmed the results 
of Refs. H[ and were able to calculate the binding energies for the three-boson system to 
extremely high accuracy. 

In this paper, we consider the four-body system with short-range interactions. The four- 
body problem has previously been studied in variety of approaches. Early studies include 
the Yakubovky equations for local potentials using the Hilbert-Schmidt expansion [14] , the 
Schrodinger equation with separable two-body potentials [l5j], and field-theoretical models 
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with separable expansions of the three-body T-matrix The four-body problem of 4 He 
atoms was investigated by Nakaichi-Maeda and Lim using the Yakubovsky equations with 



a unitary pole approximation for the S-wave (2) + (2) and (3) + (T) subamplitudes |17[. For 
a review of these and other early studies see, e.g., Refs. jig, Il9f. For an overview of recent 
calculations for the four-body system of 4 He atoms, see Refs . |2ll I22I ] . A general review of 
theoretical studies of small 4 He clusters can be found in Ref. [23j. In Ref. |2J], a benchmark 
calculation comparing various modern calculational approaches to the nuclear four-body 
problem was carried out. 

The purpose of this paper is to study the four-boson system with short-range interactions 
and large scattering length in an effective theory. We will work at leading order in l/\a\ and 
use the framework of non-relativistic quantum mechanics to construct an effective interaction 
potential. This approach has the advantage that one can immediately start from the well- 
known Yakubovsky equations for the four-body system j25[. The four-boson binding energies 
are obtained by solving the Yakubovsky equations for the effective interaction potential. The 
solution of the four-boson problem in effective theory is important in several respects: 

First, it can immediately be applied to the atomic problem of 4 He atoms and is a first 
step towards the four-body problem in nuclear physics which is complicated by spin and 
isospin. The scattering length of 4 He atoms a m 100 A is much larger than its effective 
range r e « 7 A which can be taken as an estimate of the natural low-energy length scale 
/. 4 He atoms are therefore an ideal application for our theory. The three-body system of 



4 He atoms has been investigated in Refs. ja |26|, |2jj using effective field theory. While the 
universal properties of the three-body system of 4 He atoms were discussed in Ref. (26| , this 
has not been done for the four-body system. 

Second, the renormalization of the four-body system in an effective theory is an open 
question. It is clear that low-energy four-body observables must depend on a two-body 
parameter and a three-body parameter. However, it is not known whether a four-body 
parameter is also required to calculate low-energy four-body observables up to corrections 
suppressed by l/\a\. 

The theoretical situation concerning this question appears confusing. On the one hand, 
there is a renormalization argument for 5-function pair potentials that indicates that a 
new four-body parameter is required to calculate four-body binding energies |2a| . On the 
other hand, Amado and Greenwood have evaluated the trace of the four-body kernel and 
concluded that the Efimov effect is absent in the four-body system j2^|. This result suggests 
that a four-body parameter should not be necessary at leading order in l/\a\. There is 
some circumstantial evidence in favor of the latter possibility from the four-body problem in 
nuclear physics. There is a correlation called the "Tjon line" between the binding energy B t 
of the triton and the binding energy B a of the a particle 0|. Calculations of these binding 
energies using modern phenomenological nucleon-nucleon interaction potentials give results 
that underestimate both binding energies but cluster along a line in the B t -B a plane. By 
adding a three-body potential whose strength is adjusted to get the correct value for B t , 
one also gets an accurate result for B a (cf. Ref. [30]). This conclusion also holds for chiral 
nuclear potentials derived from an effective field theory with explicit pions The aim of 
the present work is to study the question about the requirement of the four-body interaction 
by an explicit calculation in the controlled environment of an effective theory. 

The paper is organized as follows. In Sec. m we will review the description of the two- 
and three-body bound state problem in effective theory and extend this framework to the 
four-body bound state problem. In Sec. Illll we will discuss the renormalization and present 
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numerical results for the case of 4 He atoms. Finally, we close with a summary and outlook 
in Sec. EH 



II. FEW-BODY BOUND STATE EQUATIONS IN EFFECTIVE THEORY 

The effective low-energy interaction potential generated by a non-relativistic EFT with 
short-range interactions can be written down in a momentum expansion. In the two-body 
S-wave sector, it takes the general form 

(k'|V|k) = A 2 + A 2 , 2 (A; 2 + fc' 2 )/2 + ... , (1) 

where k and k' are the relative three-momenta of the incoming and outgoing particles, 
respectively. Because of Galilean invariance, the interaction can only depend on the relative 
momenta. Similar expressions can be derived for three- and higher-body interactions. The 
exact form of the potential depends on the specific regularization scheme used. The low- 
energy observables, however, are independent of the regularization scheme (up to higher 
order corrections) and one can choose a convenient scheme for practical calculations. 1 

In a momentum cutoff scheme, the potential in Eq. ([TJ can be regularized by multiplying 
with a Gaussian regulator function, exp[— (k 2 + k' 2 )/A 2 ], with the cutoff parameter A. This 
factor strongly suppresses high-momentum modes in the region k, k' > A where the effective 
potential is not valid. The cutoff dependence of the coefficients A 2 (A), A 2j2 (A), ... is de- 
termined by the requirement that low-energy observables are independent of A. Of course, 
the expansion in Eq. is only useful in conjunction with a power counting scheme that 
determines the relative importance of the various terms at low energy. In the case of large 
scattering length a, the leading order is given by the A 2 term which must be iterated to 
all orders, while the other terms give rise to higher-order corrections that can be included 
perturbatively p, |6(. In this paper, we will work to leading order in the large scattering 
length a and include only the A 2 term. 

In the three-body system, a momentum-independent three-body interaction term A3 must 
be included together with A 2 already at leading order 0, Hf. Without this three-body 
interaction low-energy observables show a strong cutoff dependence and the system can not 
be renormalized. Effective range effects and other higher-order corrections can be included 
as well BUI 111. 



The power counting for the four-body system has not been formulated yet. In order to 
see whether the minimal set of interactions A 2 and A3 is sufficient, we will calculate the 
bound states in the four-body system and study their cutoff dependence. We will use a 
momentum cutoff regularization scheme as described above. A strong cutoff dependence of 
the binding energies would indicate that a four-body interaction term is required. If the 
four-body binding energies are stable under variations of the cutoff, this would suggest that 
the four-body interaction is a sub leading effect. 

In order to set up our conventions and formalism, we will first review the bound state 
equations for the two- and three-body system and then set up the four-body equations. 



1 For a comparison of different regularization schemes in the nuclear two-body problem and chiral pertur- 
bation theory, see Refs. and 0, 

respectively. 
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A. The Two-Body Sector 

We write the leading order two-body effective potential in momentum space as: 

<P|V|q> = (p|<?>A 2 (<?|q> , (2) 

where A2 denotes the two-body coupling constant and q (p) are the relative three-momenta 
in the incoming (outgoing) channel. The regulator functions 

(p\g) = g(p) = exp(-p 2 /A 2 ) , (3) 

suppress the contribution from high momentum states. In the few-body literature, they are 
often called "form factors" . 

Our normalization for plane wave and spherical wave states is 

(p|p) = 5 (3) (P - P) , (P I m\p' I' m') = 5ijP 8 w 8 mm > , (4) 
and the S-wave projection of the plane wave state with momentum p is 

(poo| P 3 s (p|pO = ( J ) ,/ ( y /) . (5) 

For convenience, we will work in units where the mass M of the bosons and Planck's constant 
H are set to unity: M = h = 1. 

The interaction (J2J) is separable and the Lippmann-Schwinger equation for the two-body 
problem can be solved analytically. The two-body t-matrix can be written as |37l |: 

t{E) = \g)r(E){g\ , (6) 

where E denotes the total energy. The two-body propagator r(E) is then given by 



r(E) 



1A.-4* / d ^ 2 |^ 



(7) 



A two-body bound-state appears as a simple pole in the two-body propagator r at energy 
E = —B 2 . Thus the two-body coupling constant \2(B 2 , A) can be fixed from the two-body 
binding energy B 2 , which is directly related to the scattering length by a = l/y/B~2 at 
leading order in l/\a\. The integrals appearing in the propagator can be expressed through 
the complementary error function 

erfc(x) = 1 - _ / e ~ t2 dt , (8) 
V n Jo 

and for E < 0, we obtain for the inverse propagator: 



t(E)- 1 = 2tt 2 



^ exp (^r) erfc ( : nr) ~ ^~^ exp (~kP) erfc (~ jt~) 



• (9) 



A similar expression for the propagator can be obtained for positive energies by adding a 
small imaginary part to E, but will not be required for our purposes. 
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B. The Three-Body Sector 



The low-energy properties of the three-body system for a given effective potential can 
be obtained by solving the Faddeev equations [38[. Faddeev's idea was to decompose the 
full three-body wave function \I/ into so-called "Faddeev components" in order to avoid 
the problem of disconnected contributions in the three-body scattering problem. For the 
three-body problem with two- and three-body interactions, the full wave function can be 
decomposed into four components |39j|: one for each two-body subcluster and one for the 
three-body cluster. 2 For identical bosons, the three-body wave function is fully symmetric 
under exchange of particles and the Faddeev equations simplify considerably. In this case, 
one only needs to solve equations involving one of the two-body Faddeev components and 
the three-body component. The two remaining two-body components can be obtained by 
permutations of particles. For more details on the Faddeev equations, we refer the reader 
to the literature |37l . liof . 

We follow Glockle and Meier (3^ and decompose the full three-body wave function as 



where P = P 13 P 23 + P 12 P 23 



(10) 



is a permutation operator that generates the two not explicitly included Faddeev components 
from i/j. The operator Pjj simply permutes particles i and j. The Faddeev equations for ip 
and ip 3 in operator form are then 



i) = G tP^ + G t^ 3 
T/>3 = G t 3 (l + P)^, 



(11) 



where Go denotes the free three-particle propagator, t is the two-body t-matrix for the two- 
body subsystem described by the component ip. t 3 is the three-body t-matrix defined by the 
solution of the three-body Lippmann-Schwinger equation with the leading order three-body 
effective interaction 

V 3 = |0A 3 (f I , (12) 
only. Since V 3 is separable, we can solve for t 3 exactly and obtain 



t 3 (E) = |£>t 3 (£)<£| , where r 3 {E) 



^3 



(13) 



The three-body regulator function |£) will be specified later. Note that t 3 is only a technical 
construct that is generally cutoff dependent and not observable. The physical three-body 
t-matrix always includes both two- and three-body forces. Since we are interested only in 
the binding energies and not in the wave functions, we can eliminate the component ip 3 and 
obtain 

ip = G tPi) + G tG t 3 (l + P)i/; . (14) 

The component ip 3 can easily be recovered by using the second line of Eq. (jll|) . 

We now derive an explicit representation of Eq. (|14|) in momentum space. We will il- 
lustrate this procedure by showing some details for the first term on the right-hand side 
of Eq. (|14|). The extension to include the second term containing three-body interaction is 



2 Note, however, that other decompositions involving only three Faddeev components are possible as well. 
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straightforward and we will only quote the final result. Furthermore, in order to understand 
the renormalization of the three-body problem it is instructive to consider the case without 
a three-body force first. The natural Jacobi momenta are given by: 



ui = -(ki - k 2 ) 



u 2 



1 



(ki + k 2 



(15) 



Here and in the following sections, we will only take S'-waves into account. Hence, we can 
project all operators accordingly and define the S'-wave projection operator: 



dUi u\ du 2 u\ \UiU 2 ) {u\U 2 \ = / T>U U 2 \U1U2) {U\U2\ ■ 



(16) 



Using the definition (uiw 2 |-0) = ip(ui, u 2 ), we can write the Faddeev equation in momentum 
space as: 



■0(ui,w 2 ) = (uiU2\GotP\ip) = 4nGo(ui,U2)g(ui) J Vu'u' 2 J Vu"u" z 'g(u' 1 )T(E 

5{u 2 - u' 2 ) 



H 2 ) 



^VI^i'V^i'V'' 



u 



(17) 



where the factor of 47r arises from our normalization of the S-wave projected two-body 
t-matrix. The free propagator for three particles in their center of mass is given by 



GoK u 2 ) = [E — E km y x =[E-u\- \u 2 2 \ 



3„,21 _1 



(18) 



where Ej, in denotes the kinetic energy. Furthermore, it should be noted that the two-body 
propagator r in Eq. (jT7j) is evaluated at the energy in the corresponding two-body subsystem: 
E — ju' 2 2 . The permutation operator P can be written as: 



{uiU2\P\u' 1 u 2 ) = J dx 



5(U1 - 7t(« 2 , W 2 ))5('«i - 7r(w 2 , U 2 )) 

tt(u 2 , u'2) 2 ir(u' 2 , u 2 ) 2 



n(u2,u' 2 ) = \J \u\ + u' 2 2 + u 2 u' 2 x. 
Using this representation of P, we can write the integral equation as: 

/■OO /"l 

^(mi,u 2 ) = 47rGo(«i,M 2 )5 , (Mi) t(E - \u\) J du' 2 u' 2 J dx g(7r(u 2 ,u' 2 )) 
x^(7r(«2,«2),«a) • 



(19) 



(20) 



(21) 



This is an homogeneous integral equation in two variables. It can be further simplified by 
defining a new function -F(w 2 ) of only one variable via 



ip(ui, u 2 ) = G Q (ui, u 2 ) g(ui) t(E - jul)F(u 2 ) 



(22) 



leading to the integral equation 



f°° 2 f 1 

F(u 2 ) = An du 2 u 2 \ dx g 



(tt(u 2 , u 2 ))G (tc(u 2 , u 2 ), u 2 )g(Tc(u' 2 , u 2 )) 



xr(E-lu' 2 2 )F(u' 2 ) 



(23) 
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FIG. 1: The shallowest three-body binding energies indicated by the solid, dashed, and dash- 
dotted lines as a function of the momentum cutoff A. The vertical dashed line indicates the cutoff 
range in which the three-body system has exactly two bound states. The horizontal solid line 
shows the energy at which the shallowest three-body state is fixed. 



The three-body binding energies are given by those values of E for which Eq. (J2%j) has a 
nontrivial solution. By expressing the two-body coupling constant A2 in terms of the binding 
energy of the shallow two-body bound state in Eq. ©, we have already renormalized the 
two-body problem. 

The three-body system is stabilized against the Thomas collapse by the presence of the 
momentum cutoff A. No three-body force is required for this purpose. After the cutoff is 
introduced, there are no bound states with binding energies B 3 > A 2 . 3 However, the three- 
body binding energies depend strongly on the value of the cutoff A. This is illustrated in 
Fig. [T] where the shallowest three-body binding energies are indicated by the solid, dashed, 
and dash-dotted lines. The three-body system has exactly two bound states in the cutoff 
range indicated by the two vertical dashed lines. 

We now include the second term on the right-hand side of Eq. (|14j) which contains the 
three-body force. The derivation of an explicit representation in momentum space proceeds 



3 This momentum cutoff A can also be thought of as introducing a finite range r ~ 1/ A for the pair potential 
in position space. The Thomas collapse is then stopped when the bound state size is of the same order 
as the range of the potential. 



S 



as for the first term. The full equation including both terms then reads 

F(u 2 ) = 4tt J du' 2 u' 2 2 J dx g(n(u 2 ,u 2 ))G (n(u 2 ,u 2 ),u 2 )g(ii(u 2 ,U2)) 

xt(E-1u' 2 2 )F(u' 2 ) 

POO 

+ (4tt) 3 / du[u[ 2 g{u[)G Q {u^u 2 )T Z {E)i{u' x ,u 2 ) 
Jo 



X 



du'(uf I duluf £«, u' 2 ! )g{u'[)G {u'l <)r(E - f 



+ / d« / d« z / dxe(7rK J <),<)^(7r« ) <)) 



-i 



xG (7r(<y 2 '),<M£ - K ,2 )F(n 2 ' 



(24) 



where £(u!,u 2 ) = (uiU 2 |£) is defined as 



f(ui,u 2 ) = exp 



u{ + \u\ 



A 2 



(25) 



Note that the term in the exponent is the kinetic energy of the three-body system. Thus, it 
is exactly the kinetic energy flowing through the three-body interaction which is limited by 
the cutoff parameter A. This choice of the cutoff function satisfies Bose symmetry explicitly. 
The factors of 47r and (4ir) 3 arise from our normalization of the S-wave projection of t and 

The value of the three-body force is determined by the renormalization condition that 
the shallowest bound state energy is fixed as the cutoff is varied. This value is denoted by 
the horizontal solid line in Fig. ^ Depending on the value of the cutoff, the three-body force 
then must provide additional attraction or repulsion in order to keep the shallowest bound 
state energy fixed as the cutoff is varied. Thus, one three-body datum (in our case the 
shallowest bound state energy) is required as input while all other low-energy three-body 
observables can be predicted. Once the shallowest bound state is fixed, the binding energies 
of the deeper bound states will also be cutoff independent. When the cutoff is increased 
and a new bound state appears at threshold, the three-body force must turn from strongly 
repulsive to strongly attractive to satisfy the renormalization condition for the shallowest 
bound state. The additional state is then added as a deep state rather than at threshold. 
Low-energy three-body observables are not affected by the additional deep bound states. As 
a consequence, the cutoff can be made arbitrarily large in the three-body system. 

The renormalization procedure determines three-body coupling constant A3 (S3, A) 
uniquely. It was used in Refs. 0, H[ to renormalize the three-body equation derived from 
nonrelativistic effective field theory with an auxiliary field for the interacting two-particle 
state. See Ref. [28| for an earlier discussion of this renormalization method. While the part 
of Eq. (|2*4"j) resulting from two-body interactions only is very similar to the corresponding 
part of the field-theoretical equation, the part containing the three-body force (in form of 
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FIG. 2: The three-body coupling constant A3 as a function of the cutoff parameter A. The solid 
blue line shows a fit of Eq. (|26|) with c = 0.076 and L3 = 23.3\/S2 to the points for A > 245^/^2- 



r 3 ) is more complicated. The simplicity of the field-theoretical equation is due the specific 
form of the three-body interaction using an auxiliary field in Refs. 0, H| 4 

We now explicitly verify the renormalization of the three-body system as described above. 
By varying the cutoff parameter A and tuning the three-body coupling A3 such that the 
shallowest three-body binding energy stays constant, we can determine the renormalization 
group evolution of A3 numerically. For the dimensionless coupling constant A3A 4 , we confirm 
the results of Refs. P, IS Ell 13||: A 4 A 3 (A) shows a limit cycle behavior and is single-valued. 5 
This limit cycle is shown in Fig. El For large values of the cutoff A, the dimensionless three- 
body coupling constant A 4 A 3 flows towards an ultraviolet limit cycle. For A — > 00, it has 
the limiting behavior 



c sin(s ln(A/L 3 
A 3(AJ - TJ 



arctan(l/so)) 



A 4 sin(so ln(A/L 3 ) + arctan(l/so)) 



(26) 



where sq ~ 1.00624 is a transcendental number that determines the period of the limit cycle. 
If the cutoff A is multiplied by a factor exp(mr / sq) ~ (22.7) n with n an integer, the three- 
body coupling A3 is unchanged. L3 is a three-body parameter generated by dimensional 



4 One would expect that using auxiliary fields would also simplify the four-body equations. However, it is 
not obvious how to treat the (2)+(2) clusters in the intermediate state without introducing uncontrolled 
approximations . 

5 Note that A3 is defined with the opposite sign of the three-body coupling constant H in the field theory 
formulation tZllg. 
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transmutation. One can either specify the dimensionless coupling constant A 4 A3 (A) and 
the cutoff A or the dimensionful three-body parameter L 3 . The constant c is universal and 
independent of L 3 within our numerical accuracy. We have determined c numerically by 
fitting A 4 A 3 (A) for different three-body parameters L 3 to Eq.()26|) and found: 

c = 0.074 ±0.003, (27) 

where the error has been estimated from the observed variation in the fit results for c. The 
three-body parameter L 3 can be determined by fixing a three-body binding energy B 3 . Of 
course, one could also use a three-body binding energy directly to characterize the value of 
the three-body coupling A 3 at a given cutoff. However, it is advantageous to use L 3 because 
Eq. (|26|) takes a particularly simple form in terms of L 3 . 

One might expect that the smooth Gaussian regulator function we use would lead to 
better behaved numerical solutions of the three-body equations than the sharp cutoff used 
for the field-theoretical equation in Refs. 0, H| . However, it turns out that the convergence 
of the three-body coupling to the ultraviolet limit cycle is significantly slower than in the 
effective field theory formulation. This is due to the more complicated structure of the part 
of Eq. (}2"4*j) containing the three-body force. For very large loop momenta, the three-body 
part of Eq. ()24|) simplifies considerably and has a similar form as in the EFT formulation. 
In this limit , it is possible to derive the general form of Eq. (J26j) . 

In g eneral, we recover the results for the three-body binding energies from Refs. 0, HI 
GIElJ. However, we note that somewhat higher cutoffs are required to reach convergent 
results for the deeper three-body bound states. The numerical effort for the solution of the 
three-body equations becomes larger as the value of the three-body coupling A3 is increased. 
This is also related to the relatively complicated structure of the part of Eq. (j24j) containing 
the three-body force. 



C. The Four-Body Sector 

We now turn to the four-body sector. The four-body binding energies are given by the 
non-trivial solutions of the Yakubovsky equations |25[ which are based on a generalization of 
the decomposition into Faddeev components for the three-body system. The full four-body 
wave function \1/ is first decomposed into Faddeev components, followed by a second decom- 
position into so-called "Yakubovsky components" . In the case of identical bosons, one ends 
up with two Yakubovsky components ipA and ips- We start from the Yakubovsky equa- 
tions including a general three-body force in the form written down by Glockle and Kamada 



41[. The full four-body bound state wave function is decomposed into the Yakubovsky 



components if) a and i^b via 

* = (1 + (1 + P)P 34 )(1 + P)^a + (1 + P)(l + P)^B , (28) 
where exchanges particles % and j, P is defined in Eq. (II Oj) , and P is given by 

P = P 13 P 24 • (29) 
The equations for the two wave function components read: 

ipA = G t 12 P[(l + P U )^ A + M + ~(1 + G t 12 )G V 3 V > 

= G t 12 P[(l + P 34 )ip A + ips] , (30) 



II 



where t± 2 denotes the two-body t-matrix for particles 1 and 2 and V 3 is the three-body 
force defined in Eq. (|12j). Note that the three-body force couples to the full four-body wave 
function The factor of one third in front of the three-body force term arises because we 
insert the full three-body interaction for V 3 . This is possible since we consider three-body 
contact interactions which are symmetric under the exchange of any pair of particles. 

In order to describe the four-body system at rest, three Jacobi momenta are required. 
The structure of the four-body equations is more complex than that of the three-body 
equations because both (3) + (l) and (2) + (2) fragmentations can occur. As a consequence, 
two different sets of Jacobi momenta are required. The (3) + (l) fragmentation is described 
by the vector 

u 3 = ^[k 4 - ^(kx + ka + ka)] , (31) 

along with the Jacobi momenta for the three-body system given in Eq. ()15j) . The (2) + (2) 
fragmentation is described by the set: 



Vi 



2 V 1 



k 2 ) , v 2 



2( k i + k 2) 



2^3 + k 4 ) 



and 



v 3 



k 4 ) . (32) 



Depending on which Yakubovsky component and operator in Eq. (J30|) is concerned, one of 
the two sets of coordinates will be more convenient than the other. 

As in the previous subsection, we will only display the derivation of the analytic expression 
for the Yakubovsky equations without the three-body force. It is natural to evaluate the 
Yakubovsky component i[) A as a function of Ui,u 2 , and u 3 : 



(uiU 2 U 3 \lp A ) = 4<7T G (ui,u 2 ,u 3 )g(ui)r(E 



4 U 2 3 "3; 



X 



du' 2 u 2 



/2 



dx g{Tr(u 2 , u 2 )){ir(u' 2 , u 2 )u 2 u 3 \ip A ) 



+ - / du 2 u' 2 dx dx' g(n(u2,u 2 )) 



x (tt(Y 2 , u 2 )u 2 (u 2 , u 3 )u 3 (u 2 , u 3 ) \ip A ) 



+ - 



-1 



du 2 u 2 dx dx' g(Tc(u2,u 2 )) 



-1 



x (ir(u 2 , U2)v 2 (u' 2 , u 3 )v 3 (u' 2 , u 3 ) \ip B ) 



(33) 



where the two-body propagator t(E — ju 2 — |m|) is again evaluated at the energy of the 
corresponding two-body subsystem. G (ui, u 2 , u 3 ) is the free four-particle propagator in the 
center-of-mass system of the four particles and tt(u2,u' 2 ) is defined in Eq. (j^U)). The second 
Yakubovsky component ipB is computed as a function of the momenta Vi,v 2 , and v 3 : 



(viv 2 v 3 \4> B ) = 4vr G (v 1 ,v 2 , v 3 )g(v 1 )r(E - \v\ - v 3 ) 



d^ 3 v' 3 2 g(v' 3 )(v 3 V2v' 3 \^ B ) 



- I dv 3 V 3 I dx 9(V' 3 )(V 3 U2{V2,V' 3 ) U 3 (v 2 ,V 3 )\lp A ) 



(34) 
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In Eqs. (jHm El|) ? we have used the abbreviations 

U 2 {u 2 ,U 3 ) = a/ \ u l + ff U l + \l U ^ X ' ' U 3 {U 2 ,U 3 ) 



Ui 2 

9 "3 



\U 2 U 3 X 



V 2 (U 2 ,U 3 ) = \Ju 2 2 + |tt§ + |m2M3X' , f 3 (w 2 ,W 3 ) = ^ |«1 + |tt§ ~ §M 2 M3^' 



^2(^2,^3) = yj \v\ + §U§ - §f 2 ^' , M 3 (W2,^3) = 1/ 



(35) 



Similar to the three-body case, Eqs. (|H3l EMJ) can be simplified from two coupled integral 
equations in three variables to two coupled integral equations in two variables. Performing 
the substitutions 



(uiu 2 u 3 \ij A ) = g(u 1 )G (u 1 ,u 2 ,u 3 )F A (u 2 ,u 3 ) 
(viv 2 v 3 \i/)b) = g(v 1 )G (v 1 ,v 2 ,v 3 )F B (v 2 ,v 3 ) , 

we obtain the integral equations for F A (u 2 , u 3 ) and F B (v 2 , v 3 ): 

F a (u 2 ,u 3 ) = 4vr t(E - \u\ - §u§) 



du' 2 u' 2 / dxg(ir(u 2 , u 2 ))g(7r(u 2 ,u 2 )) 



(36) 



xG (n(u 2 , u 2 ), u 2 , u 3 )F A (u' 2 , u 3 ) 



+ - J du' 2 u 2 J dx J dx'g(7r(u 2 ,u 2 ))g(ir(u 2 ,u 2 )) 

xG Q (ir(u' 2 , u 2 ),u 2 (u 2 , u 3 ),u 3 (u 2 , u 3 ))F a (u 2 (u 2 , u 3 ),u 3 (u 2 , u 3 )) 

+ 2 J du 2 u ' 2 2 J dx J dx' g(n(u 2 ,u 2 ))g(n(u 2 ,u 2 )) 

xGo(vt(m 2 , u 2 ),v 2 (u' 2 , u 3 ),v 3 (u 2 , u 3 ))F b (v 2 (u' 2 , u 3 ),v 3 (u 2 , u 3 )) 



,(37) 



F B (v 2 ,v 3 ) = A<kt{E-\v 2 -v 2 ) 



dv' 3 v 3 2 g(v' 3 )g(v 3 )G (v 3 , v 2 , v' 3 )F B (v 2 , v' 3 ) 



- I dv' 3 v' 3 I dx g(v' 3 )g(v 3 )G (v 3 ,u 2 (v 2 ,v' 3 ),u 3 (v 2 ,v 3 )) 

xF A (u 2 (v 2 ,v' 3 ),u 3 (v 2 ,v' 3 )) 



(38) 



The inclusion of the three-body force term is straightforward but lengthy and the corre- 
sponding expressions are given in the appendix. 

In order to obtain the four-body binding energies, we have to solve the Yakubovsky 
equations with the three-body force term. The binding energies can be found by discretizing 
the above equations and calculating the eigenvalues of the resulting matrix. They are given 
by the energies at which any eigenvalue of the matrix is equal to one. The wave function is 
then given by the corresponding eigenvector. 



13 



The renormalization analysis of the four-body system is complicated by the cutoff depen- 
dence of the number of bound states in the three-body subsystems. The further the cutoff A 
is increased, the more three-body bound states appear. While the spurious deep three-body 
states have no influence on low-energy three-body observables, they create an instability in 
the four-body system which can collapse into a deep three-body bound state plus another 
particle. This limits cutoff variations to an interval A < A < 22.7 A for some A , in which 
the number of three-body bound states remains constant. Since the cutoff can still be varied 
by more than a factor of ten, we are nevertheless able to study the renormalization proper- 
ties and obtain converged numerical results. Alternatively, one could explicitly subtract out 
the spurious bound states from the three-body t-matrix. We will come back to this question 
in the next section. 



III. RENORMALIZATION AND NUMERICAL RESULTS 

In this section, we will discuss the renormalization of the four-boson system and present 
some numerical results for the four-body system of 4 He atoms. For convenience, we will set 
Boltzmann's constant to unity: k = 1. Since the scattering length of 4 He atoms is much 
larger than their effective range, they are an ideal application for our theory and a leading 
order calculation should be accurate to about 10% (since I /a ~ 10%). 

The quantitative experimental information on low-energy 4 He atoms, however, is rather 
limited. Using diffraction of a molecular beam of small 4 He clusters from a transmission 
grating, the bond length of the 4 He dimer has been measured to be (r) = (52 ± 4) A j4^| ■ 
This value is an order magnitude larger than their effective range r e « 7 A, which can 
be taken as an estimate of the natural low-energy length scale /. The scattering length 
a = (104^ 8 )A and the dimer binding energy B 2 = (l.ltnj} mK were derived from the 
measured bond length using the zero range approximation [42j. The 4 He trimer, tetramer, 
and several larger 4 He clusters have been observed 4^, 44]| , but no quantitative experimental 



information about their binding energies is available to date. 

However, there is a large number of theoretical calculations using realistic 4 He potentials 
for the trimer ( 4 He3). These calculations typically predict a trimer ground state with an 
energy of about 120 mK and one excited state with a binding energy of about 2 mK 0, 



46l. 1471 l48j| . The ground and excited states of the tetramer ( Hezi) and larger clusters have 
been calculated by Blume and Greene (BG) jij • They have used the LM2M2 potential ji^ 
and a combination of Monte Carlo methods and the adiabatic hyperspherical approximation 
Their results for the trimer energies agree with the exact three-body calculations of Refs. |4~ 

HESS. 



In the absence of quantitative experimental information on the three-body clusters, we 
take the binding energy of the 4 He trimer excited state from theoretical calculations using 
the LM2M2 potential as input to fix L 3 . We use the value of BG: B^/B 2 = 1.767 for this 
purpose (2ol . Elf . We can then calculate the trimer ground state and the tetramer binding 
energies for the LM2M2 potential based on low-energy universality. Before we present our 
numerical results for the 4 He4 system, we discuss the renormalization of the four-body 
problem. 

In Fig. |2J we have plotted the three-body ground state energy B% and the four-body 
energies as a function of the cutoff A. As in the case of the trimer, the 4 He tetramer has a 
ground state B± and one excited state B± . The cutoff dependence of B^ must flatten out 
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FIG. 3: Binding energies of the three- and four-body system as a function of the cutoff A. B„, 
and Bn^ denote the ground and first excited state of the n-body system. 



(0) 



and reach a plateau as A is increased since the three-body system was renormalized by fixing 
the trimer excited state binding energy (cf. subsection IIIB|) . However, Fig. El shows 
that the four-body binding energies B^ and B± also reach a plateau as A is increased. 
The excited state energy has a negligible cutoff dependence already at fairly small 
cutoffs. For the ground state Bf\ the situation is somewhat more complicated and the 
cutoff dependence of Bf' reaches a plateau only at the largest cutoff values calculated. The 
residual cutoff dependence is about 2% for the excited state and 5% for the ground state. 
A slower convergence for the ground state is expected since the value for B^ is a factor 
four larger than for the excited state and finite cutoff effects of the order y^\E\/A are more 
important. The residual cutoff dependence for both states is at least a factor 2 smaller 
than the corrections from higher orders in the expansion in I /a which are expected to be of 
the order l/a ~ 10%. We speculate that higher precision could be achieved by increasing 
the cutoff further. As already noted earlier, this would create unphysical three-body bound 
states which create an instability in the four-body system. These states would have to be 
subtracted explicitly. While such a subtraction is possible, this is beyond the scope of our 
paper and we will not attempt such a subtraction here. 
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system 


[mK] 


BW [mK] 


4g t mK i 


4g i mK i 


4 He 3 


127 


[2.186] 


125.5 


2.186 


4 He 4 


492 


128 


559.7 


132.7 



TABLE I: Binding energies of the 4 He trimer and tetramer in mK. The two right columns show 
the results by Blume and Greene |2ll | (denoted by the index BG) while the two left columns show 
our results. The number in brackets was used as input to fix L3. 

Taken together, the above observations provide strong numerical evidence that the four- 
body binding energies are cutoff independent up to higher order corrections in I /a. In 
particular, a four-body force with limit cycle behavior would lead to a much stronger cutoff 
dependence of the binding energies (cf. Fig. ^) and can be excluded. The occurrence of 
the plateaus for Bf^ and B± in Fig. |3] suggests that a four-body force is not required 
for renormalization of the four-body system at leading order in I /a. Renormalization of 
the three-body system automatically generates cutoff-independent results for the four-body 
binding energies. As a consequence, the four-body binding energies can be predicted from 
two- and three-body input alone. 

We now turn to our numerical results for the four-body system of 4 He atoms. From the 
plateaus in Fig. El we can read off the values of the binding energies. A comparison of 



our results with the values obtained by BG 2l| is shown in Table H] The results of their 



calculation for the trimer and tetramer are given in the two right columns of Table HI while 
our results are given in the two left columns. In general, our results are in good agreement 
with the values of BG. For the trimer ground state and the tetramer excited state energies, 
we obtain the values B$ = 127 mK and B± = 128 mK, respectively. For the tetramer 
ground state, we obtain B^ = 492 mK. While the value of is already relatively large, 
it is still a factor three smaller than the natural four-body energy scale ~ 1.5 K where the 
effective theory description is expected to break down. 

The natural energy scales can be estimated as follows: For two particles, it is directly 
determined by the natural length scale I and the mass M of the particles: 62 ~ ti 2 /(Ml 2 ) rs 
250 mK. For three and four particles, this estimate should be scaled according to the number 
of pairs available, 6 leading to the values ~ 750 mK and e 4 ~ 1.5 K. This estimate can 
be made for cutoff values at which the three-body force vanishes. Since all observables are 
independent of the cutoff, however, it is valid for arbitrary cutoffs. Our short-range effective 
theory can describe n-body bound states with binding energies B n <C e n . For deeper bound 
states closer to the natural energy scale the errors are expected to increase. 

The values in Table have been computed at a cutoff of A = 235 which is close 
to the largest possible value with only two three-body bound states. Our values for B^ 
and B^ agree with the BG values to within 12% and 3%, respectively. The dominant 
correction to our results is due to effective range effects which are not included in our 
leading order calculation. These deviations are within the expected accuracy of the effective 
theory. We expect the effective range corrections to the leading order result to be of the 



6 We are grateful to Eric Braaten for suggesting this scaling of the natural energy scale according to the 
number of pairs to us. 
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order r e /a ~ 10%. From the residual cutoff dependence, we estimate the numerical error of 
our calculation to be of the order of 2% for the excited state and 5% for the ground state. 
For fixed value of the cutoff A, our calculations are numerically accurate to 3 significant 
digits. 

The large scattering length of 4 He atoms also leads to universal properties in the four-body 
problem. A typical example is the existence of correlations between different observables. 
These correlations become manifest in universal scaling functions relating dimensionless 
combinations of observables. Various scaling functions for the three-body system of 4 He 
atoms were studied in Refs. EHEJ. Here, we consider the four-body binding energies as 



a function of the three-body binding energies. In phenomenological calculations of the four- 
nucleon system, an approximately linear correlation between the three- and four-nucleon 
binding energies for various nucleon-nucleon potentials is observed: the Tjon line This 
correlation is approximately linear for the relevant range of binding energies [3(J . Since the 
4 He trimer and tetramer have both a ground state and an excited state, there are four "Tjon 
lines" in this case. 

In Fig. El we show the correlations between the ground and excited state energies of the 
4 He trimer and tetramer. The two plots in the upper row show the energies of the tetramer 
excited state B± (left panel) and ground state B± (right panel) as a function of trimer 
excited state energy B^ , while the lower row shows the same quantities as a function of the 
trimer ground state energy B% . The solid line is the leading order result of our effective 
theory calculation and the cross denotes the result of the calculation by BG for the LM2M2 
potential |2l|]. For the ground states of the trimer and tetramer, calculations with other 4 He 
potentials are available as well. As an example, we show the results for the TTY, HFD-B, 
and HFDHE2 potentials taken from Refs. O,!^. 

Similar to the nuclear sector, we find an approximately linear correlation over the range of 
binding energies relevant to 4 He atoms. The calculations for the realistic 4 He potentials fall 
close to the universal scaling curves from our effective theory. For the correlation between the 
ground state energies, the "Tjon line" is directly evident in the potential model calculations 
shown in Fig. 0J If calculations with other potentials were performed for the excited state 
energies, they would also fall on a line parallel to the universal scaling curve from effective 
theory. The deviation of calculations using realistic potentials from the universal line is 
mainly governed by effective range corrections which are expected to be of the order 10% 
for 4 He atoms. For the tetramer ground state, this deviation is about a factor four larger 
than for the excited state. 

While the correlations in Fig. 0] are approximately linear over the range of binding ener- 
gies calculated, we expect them to become nonlinear for a sufficiently large range of binding 
energies. Similar nonlinearities were observed for correlations between three-body observ- 
ables in Refs. |H|26|. Our results suggest that the "Tjon lines" are universal properties of 
few-body systems with short-range interactions and large scattering length. They do not 
depend on the details of the short-distance physics which are very different in atomic and 
nuclear systems. It is interesting to note that in nuclear physics only a correlation between 
the ground-state energies of the three- and four-body systems has been observed. We ex- 
pect this correlation to hold for all S-wave states that are within the range of validity of 
an effective theory with contact interactions. It would be interesting to see whether such 
a correlation also holds for the excited J p = + state of the a particle above the p+ 3 H 
threshold. 
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FIG. 4: 



The correlations between the ground and excited state energies of the 4 He trimer and 

' ■* (left panel) and ground state energy 
(i) 

3 ' 



tetramer. Upper row: the four-body excited state energy B 4 
? (o) 



B\ u ' (right panel) as a function of the three-body excited state energy B 



Lower row: the 

same quantities as a function of three-body ground state energy B% . The solid line shows the 
leading order effective theory result and the cross denotes the calculation for the LM2M2 potential 
by Blume and Greene 2]|. The triangles show the results for the TTY, HFD-B, and HFDHE2 
potentials 



17 



We have fitted the scaling functions shown in Fig. 0] with linear expressions and obtained: 

d(0) d(0) o(0) 

—J- = -24.752 + 4.075 -J- , 69 < -5- < 142 , (39) 

B~2 B 2 B 2 

d(0) d(1) 

— = -742.0 + 645.1 -J- , 1.54 < -J- < 2.00 , (40) 

B 2 -D2 -D2 

o(l) d(0) o(0) 

= -0.662 + 1.034 -J- , 65<^-<125, (41) 

£> 2 -D2 -D2 
dW 

4 1 70 n 1 1 Kfi A 3 iko^3 



-178.0 + 159.4 , 1.52 < -±- < 1.92 . (42) 

B 2 B 2 B 2 

These relations can be used to predict the tetramer ground and excited state energies for any 
potential for which one of the trimer energies and the dimer binding energy are known. The 
expressions (j3HJ)-(J32!) are of the same accuracy as our explicit calculations (see the discussion 
above). They are expected to be most accurate for the excited states. 
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IV. SUMMARY AND OUTLOOK 



In this paper, we have studied the four-body system with short-range interactions and 
large scattering length. We have concentrated on the bound state problem of four bosons 
starting from the Yakubovsky equations j25[. We have constructed an effective interaction 
potential including both a two- and three-body contact interaction. This is the minimal set 
of contact interactions required for renormalization of the three-body problem 0, H[. The 
two parameters of the effective potential were determined from matching to the binding 
energy of the dimer and the excited state of the trimer. We have then solved the four-body 
bound state problem under the assumption that no four-body interaction is required for 
renormalization at leading order. We found that after renormalizing the two- and three- 
body subsystems, the four-body binding energies were automatically independent of the 
ultraviolet cutoff. This result suggests that the four-body interaction is not of leading order 
and the low-energy four-body observables are determined by properties of the two- and 
three-body systems up to corrections suppressed by l/\a\. Although we have considered 
only the four-boson bound state problem with large scattering length explicitly, we expect 
this result to hold for all low-energy four-body observables. 

We have applied this effective theory to the four-body system of 4 He atoms and calculated 
the ground and excited state energies of the 4 He tetramer. In the absence of experimental 
information on the 4 He trimer, we have taken the excited state energy of the 4 He trimer 
as calculated by Blume and Greene for the LM2M2 potential [2l( as input to determine 
the three-body parameter L 3 . For the binding energies of the trimer ground state and the 
tetramer ground and excited states, we find the values = 127 mK, B± = 128 mK, and 
B^ = 492 mK, respectively. The latter energy is still about a factor three smaller than the 
natural four-body energy scale ~ 1.5 K where the effective theory is expected to break down. 
Our values for and B± agree with the calculation by Blume and Greene to within 12% 
and 3%, respectively. These deviations are consistent with the expected accuracy at leading 
order in the large scattering length of about l/\a\ ~ 10%. 

The large scattering length of 4 He atoms leads to universal properties such as universal 
scaling functions. We have calculated the universal scaling functions relating the tetramer 
energies to the trimer energies. The correlations are approximately linear in the region 
of binding energies relevant for 4 He atoms. As expected from low-energy universality, the 
results of various calculations using realistic 4 He potentials fall close to the universal scaling 
curves. Corrections to the scaling curves are mainly governed by effective range effects. We 
have fitted the calculated scaling functions with linear expressions (|39|) - (j42|) that can be 
used to obtain the tetramer binding energies at leading order in l/\a\ for any potential if one 
of the trimer binding energies is known. 

There are a number of directions that should be pursued in future work. While we have 
demonstrated that a four-body force is not necessary to renormalize the four-body system 
to leading order, the general power counting for four-body forces is still not understood. At 
which order does the leading four-body interaction enter? In the three-body system, e.g., 
the first order correction is due to the two-body effective range. If a similar situation holds 
in the four-body system, it would be possible to predict low-energy four-body observables 
up to corrections of order (I /a) 2 from two- and three-body information alone. 

The extension of the effective theory to calculate four-body scattering observables would 
be very valuable. The knowledge of the dimer-dimer scattering length, for example, is 
important for experiments with ultracold atoms. For the simpler problem of fermions with 
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two spin states (where the three-body parameter L 3 does not contribute), the dimer-dimer 
scattering length was recently calculated |52T |. 

Whether this effective theory can be applied to the nuclear four-body system like the 
pionful theory jU is an open question. While it is straightforward to generalize the effective 
theory to include spin and isospin, it is not clear whether an effective theory without explicit 
pions will be adequate for the a particle ground state with a binding energy of about 28 
MeV. This question deserves further study. The effective theory might also help to shed 
some light on the renewed speculations about the existence of a shallow tetraneutron bound 
state 1 531. 
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APPENDIX A: YAKUBOVSKY EQUATIONS WITH THREE-BODY FORCE 

In this appendix, we derive the analytic expressions for the three-body force term in the 
Yakubovsky equations (j3*Uj) in momentum space. The three-body force term 

^( Ul u 2 u 3 \(l + G t 12 )G V 3 \^) , (Al) 

couples to the full four-body wave function \l/ (cf. Eq. (J2EJ)) 5 which is related to the Yakuvosky 
components ipA and ips via 

* = (1 + P34 + PP34 + P + PuP + PPuP)^a + (1 + P + P + PP)^b ■ (A2) 

In order to simplify Eq. ()A1|) . we consider a term with an arbitrary permutation operator 
X acting on ip^. Inserting the S'-wave projection operator, Eq. ([16)1. we can write 

^( Ul u 2 u 3 \(l + G t 12 )G V 3 X\^ A ) = 



y(4vr) 2 G (Ml,M2,M3) 



J du[u[ 2 du 2 u 2 2 Vu'' u' 2 £(ui, u 2 )^(u[, u 



+4tt J du' l u' 2 du';ufdu^ufVu"' U '" 2 g(u l )T(E - \u\ - §«|)<?K) 
xG (u[, u 2 , u 3 )£(«'i, u 2 )i{u'{, u 2 ^(u"u 2 ^ 3 \X\u"'u 2 u 3 ")(u"'u 2 u 3 "\^ A ) 



(A3) 



where the factors of 4tt and (47r) 2 arise from the S"-wave projection of the two-body t-matrix 
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and three-body potential, respectively. This expression can be rewritten as: 
^(«i«2U 3 |(l + G t 12 )G V 3 X\i/> A ) = 



A 



^(4tt) 2 G'o(mi, u 2 , u 3 ) C(ui, u 2 ) + 47r#(w 1 )r(£ - \u\ - |m|)X(m 2 , u 3 ) 



x / du / 1 u' 2 du / 2 u' 2 2 Vu"u" 2 ^(u' 1 , u' 2 ){u[u' 2 u 3 \X\u'[ulul)(u'[u'iul\i) A ) 



A 



3 (47t) 2 G ( Mi , u 2 , u 3 ) [t( Ul , u 2 ) + 47ig( Ul )r(E - \u\ - §«i) J(« 2 , u 3 )\ K^\u 3 ) (A4) 
where we have defined the quantities 

Z(u 1 ,u 2 ) = J du[u[ 2 G (u[,u 2 ,u 3 )g(u' 1 )£(u' 1 ,u 2 ) , (A5) 



and 



K l x\u 3 ) = J du[u[ 2 du 2 u 2 2 Vu"u" 2 Z(u[,u 2 ){u[u 2 u 3 \X^ (A6) 

An analogous expression K x can be derived for permutation operators which act on the 
second Yakubovsky component ipB- 

We proceed by giving the analytical expressions for the and which appear in 
the computation of the three-body force term. There are six combinations of permutation 
operators acting on ip A : 

Ki» = /d^durf^^yWuWW , (A7a) 

K { p^ = \J dui«i 2 dt4w2 2 ^(wi,t4) J dx' (u[u 2 (u' 2 ,u 3 )u 3 (u' 2 ,u 3 )\i! A ) (A7b) 

K pL = \j du^du'^f J\x J\x' ^n(u' 2 ,u^,u' 2 ) 

x(n(u 2 \u 2 )u 2 (u 2 \u 3 )u 3 (u 2 ,u 3 )\i) A ) , (A7c) 
K ( P A) = J du'^du'iuf y dx^TT^^O.^XTK^Kwsl^A) , (A7d) 

^Psip = ^/ du' 2 u' 2 du' 2 u' 2 J dx J dx' £(7r(u 2 (u' 2 ,u 3 ),u 2 ),u 2 ) 

x (7r(«2, m 2 (m 2 , u 3 )) u 2 u 3 (u 2 , u 3 \^ a ) , (A7e) 
^pp 34 p = ^/ du' 2 u 2 2 du 2 u 2 2 J dx J dx J dx" $,(ir(u 2 , u 2 ), u' 2 ) 

x (uj [7r(M 2 ', u' 2 ),u 2 (u 2 , u 3 )]u 2 [ir(u 2 , u' 2 ),u 2 (u 2 , u 3 )]u 3 (u 2 , u 3 ) \tp A ) , (A7f ) 
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and four combinations of operators acting on ips- 

k[ B) = \J du'^du^u'^ J dx' £(u[,u' 2 )(u[v 2 (u' 2 ,u 2 )v 3 (u' 2 ,U3)\tPb} , (A8a) 

Kj& = \j du' 2 u' 2 2 du'^uf J dx j dx' tt«<),u' 2 ) 

x(^(u^u' 2 )v 2 (u' 2 \u 3 )v 3 (u'^u 3 )\iPb) , (A8b) 
KpB) = \) du i u 'i 2du 2 u 2 2 J dx' £(u[,u 2 )(v 3 (u 2 ,u 3 )v 2 (u 2 ,u 3 )u 1 \?Pb) , (A8c) 

K { pl = \J du' 2 u' 2 du 2 u' 2 J dx J dx' £,{tt(u 2 , u 2 ), u' 2 ) 

x(v 3 (u' 2 ,u 3 )v 2 (u 2 ,u 3 )tt(u' 2 ,u' 2 )\^b) , (A8d) 
where Ui(ui,u 2 ) and W2(wi,W2) are defined as 



ui(u 1 ,u 2 ) = Jjuj + f^ul + \uiu 2 x" , m 2 (m 1; m 2 ) = Ju\ + jul - uiu 2 x" . (A9) 
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